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Abstract - We compute the decay of the velocity autocorrelation 
function, the Lyapunov exponent and the diffusion constant for the 
Sinai billiard within the framework of dynamical zeta functions. The 
asymptotic decay of the velocity autocorrelation function is found to 
be C(t) ~ c(R)/t. The Lyapunov exponent for the corresponding map 
agrees with the conjectured limit X map — > —2 \og(R) + C as R — > 
where C = 1 - 4 log 2 + 27/(2?r 2 ) • ((3). The diffusion constant of the 
associated Lorentz gas is found to be divergent D(t) ~ logi. 



1 Theory 



1.1 Chaotic averages 

Consider a chaotic Hamiltonian system with two degrees of freedom. Assign 
a weight w(xa,t) to the trajectory starting at phase space point Xq and evolving 
during time t (to point x(xo,t)) in such a way that it is multiplicative along the 
flow: w(xo,ti+t2) = w(xo, ti)w(x(x o, t%), t^)- The phase space average of w(xo, t) 
may be expanded in terms of periodic orbits as |Q 

te^'w-fepf-t ia-^a-'j/Aai ■ <" 

where r is the number of repetitions of primitive orbit p, having period T p , and 
A p is the expanding eigenvalue of the Jacobian (transverse to the flow). w p is the 
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weight integrated along with cycle p. Zeta functions are introduced by observing 
that the average (|l|) may we written as 

lim {w(x 0} t)) = lim — / e^^lds , (2) 

with the zeta function 

oo ✓ _ g rp v m+1 

p m=0 V I PI p / 

The expressions are valid also for an area preserving map, like a Poincare map 
of the above system. The time vaiable t and period T p are now discrete and usually 
denoted as n and n v . The zeta function (^) is usually expressed in the z = exp(s) 
plane. The Dirac delta in (|l|) should be replaced by a Kronecker delta and the 
integral (0) accordingly performed along a circle \z\ = C. 

We will be interested in three weights of the required type. In the discussion 
of Lyapunov exponents we will be interested in 

w x (x ,t) = A(x ,tf , (4) 

where A(xo,t) is the (leading) eigenvalue of the Jacobian of the map from xq to 
x(xo,t). It is not exactly multiplicative but we will discuss that problem below. 
We now express the Lyapunov exponent in terms of the associated zeta function 

A= lim i(log|A(xo,*)|> = lim ~(A(x ,tf} \p= 

t—KX t t^OO t dp 



1 1 f™ +e d d 
lim-—/ e st — — logZ A (s) \ p=0 ds . 
t->oo t 2m J_ iw+e dp ds 



(5) 



_im / 

In the discussion of diffusion we will consider 

W D (x ,t) = e H^o),t)-x ) ; (g) 

where x is the configuration space part of the phase space vector x. The diffusion 
constant may now be expressed in terms of the associated zeta function 

D = lim U(x(x ),t) - xo) 2 } = lim + J^ ){e ^(*oU)-*o) } | 

t^oo 2t t->oo 2t dpi dp^ 

1 1 f ioo+e st d 2 ^ 2 J 



Next assume that we have a Poincare map defined by a surface of section 
^s.o.s.- Then, using weight 

e~ s - T{xo < n) w(x ,n) , (8) 
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where T(xo,n) and w(xo,n) is the integrated time and weight along a trajectory 
starting at xq 6 ^s.o.s. an d evolving during n iterates of the map. We thus get a 
zeta function which combines the description of the map and the flow. 

oo , ^ v m+l 

2(v)=nn i-v^^ o) 

V rn=Q V l A Pl A p / 

A more lucid introduction to chaotic averages may be found in with proper 
references. 

1.2 Approximate zeta functions 

Given a time continous system, suppose we have defined a surface of section in 
such a way that the corresponding map destroys all information from one iterate 
to another. Or more precisely, the weight associated with consecutive iterates are 
uncorrelated: (wi w;+i) = (wi)(wi + i). This is called Assumption A in ref. [||. 
This can of course not be exactly fulfilled in practice. But suppose our system 
is intermittent and we put the s.o.s. in the border between the laminar and the 
chaotic phase. Then the chaotic phase will do its best to destroy the memory 
from one laminar phase to the next. We will now average the the weight over the 
surface of section. 

{e- T ^riw(x ,n)) a . o . a . = (e- sT ^w(x ,l))l o , s , , (10) 

which follows from assumption A. This type of behaviour is equivalent to their 
being only one zero zq(s) — (exp(— sT(xo, l))w(xo, l)) s .o.s. of Z{s,z) in the z— 
plane. In the language of Ruelle resonances this means that there is an infinite 
gap to the first resonance. This in turn implies that the Z(s 1 z) may be written. 

Z{s,z) = e h ^ (z-{e- sT ^^w{x ,l)) s . o . s ) , (11) 

where h(z, s) is entire in z and s. The zeta function for the flow is obtained as 
Zfiow(s) = Z(z = 1, s), see eq. (||), and we get 

Z flow {s) = 1 - (e- sT ^w(x , l)) s . . s . . (12) 

We had to put h = (z = l,s) = in order for Z(s) to obtain the correct limit 
at infinity. Normally, the best one can hope for is an finite gap to the leading 
resonance of the Poincare map. The zeta function above is then only approximate, 
indicated by putting a hat on it: Z, we also suppress the clumsy subscript flow. 
The first two terms in an expansion around s = for the special case w = 1 are 
however exact: Z(s = 0) = Z(s = 0) = and Z'(s = 0) = Z'(s = 0) for linearity 
reasons. 

We can reformulate the formula above for the special case where w(xq, 1) is a 
function of T{xq, 1) = A(xo) denoted w(A) Then the probability distribution of 
recurrence times is 

p(A) = {S(A - A(x ))) s .o. s . , (13) 
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and 

poo 

Z{s) = 1 - / w(A)p(A)e- sA dA . (14) 
Jo 

Approximate zeta functions were introduced in refs. [§], |5|, to a large extent 
inspired by ref . || , the above derivation is though much more easy to follow than 
the one in [|J. 



2 The Sinai Billiard 

We will apply the method of approximate zeta functions to the Sinai billiard[Q , 
a particle with unit speed bouncing inside a unit square with a scattering circle of 
radius R in the middle. The only sensible choice of £l s .o.s. is the disk itself. The 
disk to disk map will be hyperbolic but without a finite Markov partition, most 
likely it will have a finite gap and exhibit exponential mixing ||, In the small 
scatterer limit one may derive the following exact expression for the distribution 
of recurrence times [HOl . 



12R. 



^(2£ + C(4 - 30 log(0 + 4(£ - l) 2 log(£ - 1) - (2 - 2 log |2 - £|) £ > 1 

(15) 

where £ = A2R In the following we will e.g. need the mean < A >= l/2i? in this 
distribution, and the large A limit: p ~ l/(7r 2 i? 2 A 3 ). 



2.1 The Lyapunov exponent 

Consider a trajectory with n intersections with Q s . . s .. The consecutive disk 
to disk time of flight and scattering angles will be denoted as Ai and cti being 
functions of the s.o.s coordinate Xi. It may be shown [0 that the (by construction) 
multiplicative weight 

" 2A- 

w(x ,n) = T\— — , (16) 

- 1 1 R cos ai 

i—l 

give the correct R — ► limit of the Lyapunov exponent 

1 2A(x) 
X map = Mm -(logw(x ,n)) s . o . s + O(R) = (log- ^ T ) s . . s + 0{R) (17) 

= log(2/i?) + (log A{x)) s . . s - (logcosa(a;)) s . . s + O(R) . 

The problem is then reduced to computing the average of log A (the average of 
cos a is trivial). To that end we use eq. (15). The result is ]To| 

X map -> -2 log R + C R^oo , (18) 
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with C = 1 — 4 log 2 + 27/ (2tt 2 ) ■ C(3), confirming previous conjectures 11, 1^, [l^] 
but with C computed exactly for the first time. 

What about the Lyapunov exponent of the flow. The zeta function is now, cf. 
eqs. (|) and ©, 

Z = 1 - (A(x fe- sA ^) s . . s = s(A) + 0(s 2 logs) - (A map + 0(a)) . . . (19) 
Putting this into eq. (||) we get 

= ^l<W i+0(1 ° 6<)) = W ' <20) 

which is just Abramovs formula We got an exact result (for the weight 

([l6|)) since the relevant terms in the approximate zeta function were indeed exact. 
We saw that the leading zero also being a branch point did not do any harm to 
the Lyapunov exponent but it will in fact induce a phase transition among the 
generalized Lyapunov exponents . 

2.2 Diffusion 

Next we consider the associated Lorentz gas obtained by infinitely unfolding 
the billiard. The zeta function is, cf. eqs. (||) and ( |l2] ) 

1 r 2 " 

Z{s) = l-{e- sA e Xll3l+X2 > 32 ) = l-{e- sA — e A0cos ^d9), (21) 

2vr J 

where we used the asymptotic isotropy (which is not really necessary cf refs. |l], 
^5|), and A = \J x\ + x\. We only need to expand to second order in /3 

Z(s) = l-(e- sA )-^(A 2 e- sA )... 
= 1 -^ + 0(s 2 \og S ) + ^ ¥ (-\og S + 0(l))... . (22) 



Putting it into eq (fj) we get 



1 2 1 rSt -logs + 0(l)^_ 1 

-zoo+e 



D = ^2t^R^i * —> ^ =^(log(t) + 0(l)) , (23) 



which is indeed the exact result |L5| ]. Exact result may also be obtained for finite 
R \ ll |l5| . The exactness of the result indicates that the term — 47r ^a log s in eq. 
(p2) agrees with the corresponding one in the exact zeta function. 
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2.3 Correlation functions 



As regarding correlation functions we will restrict ourselves to observables 
A(xq) that are changed only by bounces on the disk, such as A = \v x \. The 
autocorrelation function is e.g. obtained as the time average 

C(t) = (A(ta + t)A(t ))t -(A 2 ) . (24) 

Measuring A at two points of time separated by t, then there is a probability po(t) 
that the trajectory has not hit the disk in between and according to assumption 
A we get pi 

C(t)=p (t)(A 2 ) + (l-p (t))(A) 2 -(A 2 )=p (t)V(A) , (25) 

where V(A) is the variance of A. The function po(t) may be expressed in terms of 
P (A) @ 

(26) 



The 1/A 3 decay of p(A) thus implies a 1/t of the correlation function. In fig 1 we 
compare the numerical correlation function with A = \v x \ with eq (|25|) with po(t 



computed from eq. (|l_5J). The 1/t decay law has also been suggested in refs [17 
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Figure captions 

Figure 1: Numerical correlation function (full line), for R — 0.106. The dashed 
line represents eq. (p5|) with p (t) computed from eqs. (|26]) and (|l5|). From ref . Jl(| 
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